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Abstract. I present here the main results of recent realistic, 3D, hydrodynamical 
simulations of convection at the surface of metal-poor red giant stars. I discuss the 
application of these convection simulations as time-dependent, 3D, hydrodynamical 
model atmospheres to spectral line formation calculations and abundance analyses. 
The impact of 3D models on derived elemental abundances is investigated by means of 
a differential comparison of the line strengths predicted in 3D under the assumption of 
local thermodynamic equilibrium (LTE) with the results of analogous line formation 
calculations performed with classical, ID, hydrostatic model atmospheres. The low 
surface temperatures encountered in the upper photospheric layers of 3D model 
atmospheres of very metal-poor stars cause spectral lines of neutral metals and 
molecules to appear stronger in 3D than in ID calculations. Hence, 3D elemental 
abundances derived from such lines are significantly lower than estimated by analyses 
with ID models. In particular, differential 3D— ID LTE abundances for C, N, and O 
derived from CH, NH, and OH lines are found to be in the range —0.5 to —1 dex. 
Large negative differential 3D— ID corrections to the Fe abundance are also computed 
for weak low-excitation Fe I lines. The application of metal-poor 3D models to the 
spectroscopic analysis of extremely iron-poor halo stars is discussed. 
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1. Introduction 



The first evidence of the existence of stars characterized by a lower abundance of metals 
with respect to hydrogen than in the Sun came in the Mid-20th century, in particular 
with the works of Schwarzschild & Schwarzschild (1950) and Roman (1950). These 
authors observed that giant and dwarf stars with weaker metal lines in their spectra 
typically had a larger velocity dispersion than stars with stronger spectral lines; also, 
high-velocities appeared to be associated only with the weak-line stars. These works 
suggested that population II stars (as dubbed by Baade 1944) differed from population 
I stars not only by their kinematic properties but also in terms of basic chemical 
composition, the former being characterized by an overall lower iron abundance. The 
following year Chamberlain & Aller (1951) produced the first quantitative abundance 
analysis of population II stars; the authors analysed the two stars HD 19445 and 
HD 140283 and derived significantly lower abundances of calcium and iron with respect 
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to hydrogen than in the Sun, which actually allowed to explain the peculiar spectral 
and photometric properties of these objects. 

It soon became clear that metal-poor stars were extremely interesting objects which 
could provide important clues about the formation and the chemical and dynamical 
evolution of the Galaxy. Better tools for quantitative spectroscopy of low-metallicity 
stars were therefore also required, and the first ID model stellar atmospheres of 
metal-poor stars started to appear in the late 1960s-early 1970s. Noteworthy are the 
works by Krishna Swamy (1969) and Gustafsson et al. (1975); the latter in particular 
constitutes the first extended grid of theoretical ID model atmospheres of giant stars 
with metallicities ranging from solar down to [Fe/H]= —34 Since then, ID model 
stellar atmospheres have been continuously developed and improved (see, e.g., reviews 
by Gustafsson & Jorgensen 1994, Asplund 2005), especially with regard to input physics 
and opacities, and are still nowadays the most widely used models in stellar abundance 
analyses. 

Yet, classical ID model stellar atmospheres are normally constructed under 
a number of simplifying assumptions, such as plane-parallel geometry or spherical 
symmetry, hydrostatic equilibrium, and flux constancy. Moreover, ID modelling of 
stellar atmospheres inherently relies on rudimentary implementations of convective 
energy transport, which is usually treated by means of the mixing- length theory (Bohm- 
Vitense 1958) or similar alternative formulations (e.g Canuto & Mazzitelli 1991), all 
characterized by a number of tunable but not necessarily physical free parameters. In 
late-type stars however, the convection zone reaches and appreciably influences the 
surface layers from which the emergent stellar flux is formed. For instance, high spatial 
resolution imaging of the solar photosphere immediately reveals a distinctive granulation 
pattern that reflects the bulk gas motions in the upper part of the convective zone near 
the optical surface. Given the dynamic and multi-dimensional character of convection, 
the use of time-independent, hydrostatic, ID, model atmospheres in abundance analyses 
of late-type stars is arguably a potential source of severe systematic errors. In addition, 
the strengths and exact shapes of spectral lines cannot possibly be predicted in ID 
without resorting to ad hoc fudge parameters such as micro- and macro-turbulence, 
which are supposed to account for non-thermal Doppler broadening due to bulk flows 
in the stellar atmosphere. During the past thirty years, on the other hand, various 
computer codes have been developed to perform time-dependent, 3D, hydrodynamical 
simulations of convection at the surface of late-type stars (e.g. Nordlund 1982, Nordlund 
& Dravins 1990, Stein & Nordlund 1998, Asplund et al. 1999, Freytag et al. 2002, Ludwig 
et al. 2002, Carlsson et al. 2004, Vogler 2004). Such stellar surface convection simulations 
can naturally be used as time sequences of 3D model stellar atmospheres to study 
in a self-consistent way the impact of photospheric inhomogeneities and correlated 
velocity fields on the formation of spectral lines and on spectroscopic abundance 
analyses. In fact, 3D simulations have been shown successful in reproducing several 

f [Fe/H]= log (nF e /nu) ~ log (riFe/^H)©, where np c and bh are the number densities of Fe and H, 
respectively, and the subscript symbol refers to the Sun. 
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observational constraints such as the morphology of the granulation pattern in the 
Sun as well as detailed shapes of spectral lines in other solar-like stars. Recent 
analyses based on 3D simulations of convection at the surface of the Sun, dwarfs, 
subgiants, and giants (e.g. Asplund et al. 1999, Asplund & Garcia Perez 2001, Collet 
et al. 2006, Collet et al. 2007, Asplund 2005) indicate that the structural differences 
between 3D hydrodynamical and ID hydrostatic model atmospheres can significantly 
impact the predicted strengths of spectral lines and, in turn, lead to large differences in 
the derived elemental abundances, especially at low metallicities. In the following, I will 
illustrate the main results of recent 3D simulations of convection at the surface of metal- 
poor red giants (Collet et al. 2006, Collet et al. 2007) and discuss some applications of 
these models to spectral line formation and abundance analysis. 

2. The 3D convection simulations 

The 3D, time-dependent, compressible, explicit, radiative- hydrodynamical code by 
(Stein & Nordlund 1998) has been used here to simulate convection at the surface 
of red giant stars with varying effective temperatures (T eff :=s4700 to 5100 K), surface 
gravity log g = 2.2 (cgs), and metallicity ranging from solar down to [Fe/H]= —3. 
The equations of mass, momentum, and energy conservation are solved on a discrete 
Cartesian mesh at moderate numerical resolution (100x100x125) for a representative 
volume of stellar surface large enough to cover about ten granules horizontally and 
eleven pressure scale heights vertically. In terms of optical depth in the continuum 
at A = 5000 A, the simulations extend from logT5ooo~~ 4 to logrsooo^. The upper 
and lower boundaries are open, while periodical boundary conditions are assumed 
horizontally. The simulations make use of realistic input physics: the equation of 
state comes from Mihalas et al. (1988), continuous opacities from the Uppsala opacity 
package (updated version of Gustafsson et al. 1975), line opacities from Kurucz (1992) 
and Kurucz (1993). The solar chemical composition by Grevesse & Sauval (1998) is 
adopted with the abundances of all metals scaled proportionally to the [Fe/H] value 
assumed for the simulation. The radiative heating term in the energy conservation 
equation is computed by solving the radiative transfer equation at each time-step for all 
grid-points at the surface and along eight inclined rays plus the vertical. To reduce the 
computational time, opacities are re-grouped in four opacity bins (Nordlund 1982). Also, 
local thermodynamic equilibrium (LTE) is assumed, with a Planckian source function 
and scattering treated as true absorption. 

The temperature and density structures at the surface of red giant convection 
simulations are qualitatively fairly similar to the the ones previously described in the 
works by (Asplund et al. 1999) and (Asplund & Garcia Perez 2001) on dwarfs and 
turnoff stars. Warm plasma rises from the stellar interior and rapidly cools as it 
approaches the optical surface; there it eventually becomes denser than the surrounding 
gas, turns over, and plunges back towards the interior. The bulk flows at the surface 
self-organize into a dynamical granulation pattern characterized by large, warm upflows 
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separated by a network of cool, narrow downdrafts. Figure 1 shows the spatially resolved 
outgoing intensity in the continuum opacity bin for two red giant simulation snapshots 
with similar effective temperatures, at solar and very low ([Fe/H]= —3) metallicity. 
The properties of the emergent granulation patterns are slightly different in the two 
cases: at solar metallicity, the granules are typically larger and the overall intensity 
contrast is lower than at [Fe/H]= —3. More remarkable however is the metallicity 
dependence of the temperature stratification in the 3D simulations; figure 2 illustrates, 
as a function of optical depth, the predicted temperature structures from the two 
giant simulations at [Fe/H]= and [Fe/H]= —3. The resulting thermal structures 
are compared with the ID stratifications from marcs model atmospheres (Gustafsson 
et al. 1975, Asplund et al. 1997) generated for the same stellar parameters, and with the 
same opacity data and chemical compositions as the 3D simulations. At solar metallicity, 
the mean temperature structure in the upper photosphere of the 3D simulation appears 
very similar to the stratification from the corresponding ID model where radiative 
equilibrium is enforced by the stationary condition. At very low metallicity, on the 
contrary, the temperature in the upper layers of the 3D simulations tends to remain 
significantly below the radiative equilibrium value marked by the corresponding ID 
model. The reason for this difference can be found by looking at the energy conservation 
equation. The temperature in the upper layers of 3D simulations is mostly regulated by 
two competing mechanisms: radiative heating due to the reabsorption by spectral lines 
of radiation coming from deeper in, and adiabatic cooling due to the expansion of the 
upflowing gas in granules. At very low metallicity, fewer and weaker lines contribute 
to the total opacity, therefore the significance of radiative heating from spectral lines 
is reduced, adiabatic cooling becomes preponderant, and the balance between the two 
mechanisms is achieved at lower temperatures than in stationary ID models. 

At [Fe/H]= —3, the temperature difference between the mean 3D and ID 
photospheric structures is substantial and can amount to ~1000 K. With regard to stellar 
spectroscopy, such differences, as well as temperature and density inhomogeneities, 
can have a dramatic effect on the predicted strengths of spectral lines. The cooler 
photospheric stratification of the 3D convection simulations of very metal-poor stars is 
expected to have affect significantly temperature-sensitive features, such as molecular 
lines or weak low-excitation lines from minority species (e.g. neutral metals). Gravity- 
sensitive features are also affected as the lower photospheric temperatures of the 3D 
simulations also imply reduced gas and electron pressures. 

3. Spectral line formation 

As mentioned above, 3D simulations of convection at the surface of late-type stars can 
be used as time-dependent, 3D, hydrodynamical model atmospheres for spectral line 
formation purposes. In this section, I will sketch the general strategy used to compute 
line profiles with 3D model atmospheres and to derive differential 3D— ID elemental 
abundances. From the full convection simulations, representative sequences of typically 
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Figure 1. Spatially resolved emergent intensity in the continuum bin for two snapshots 
of 3D hydrodynamical simulations of red giants at [Fe/H]= 0.0 (left) and [Fe/H]= —3.0 
(right); the characteristic surface granulation pattern is apparent. 
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Figure 2. Temperature structure of two snapshots of 3D hydrodynamical simulations 
of red giants at [Fe/H]= 0.0 (left) and [Fe/H]= -3.0 (right). Gray shaded 
area: temperature distribution as a function of optical depth in the continuum 
at A = 5000 A; darker areas indicate values with higher occurrence. Solid line: 
temperature stratification of a ID MARCS model atmosphere computed for the same 
stellar parameters. Dashed line: mean temperature stratification of the 3D snapshot 
(temperature averages taken on surfaces of constant optical depth). 



~50-100 snapshots are taken at regular intervals in time; these time-series should cover 
sufficiently long periods to provide good statistical samples of the evolution of the surface 
granulation patterns. The snapshots are interpolated to a finer depth-scale to increase 
the spatial resolution in the atmospheric layers with optical depth less than logr 50 oo~2.5 
to achieve higher numerical accuracy. Flux profiles are computed for a variety of lines 
from metals and molecules. The radiative transfer equation is solved for about 60 to 
100 wavelength-points per spectral line profile, along typically ~ 10-30 rays, and for all 
grid-points at the surface of the 3D model. As during the construction of the simulation, 
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the assumption of LTE is made, with scattering treated as true absorption. The same 
numerical code, input physics, and opacity data are used to compute spectral line profiles 
with ID model atmospheres corresponding to the same stellar parameters; in the ID 
calculations, a non-zero value of the micro-turbulence is adopted to mimic non-thermal 
Doppler broadening of the lines. I would like to emphasize once more here that in 
3D, on the contrary, no free parameters such as micro- or macro-turbulence enter the 
calculations: only the velocity fields predicted by the convection simulations are taken 
into account to reproduce non-thermal broadening, asymmetries, and wavelength shifts 
associated with the bulk motions of the gas in the stellar atmosphere. 

The impact of 3D hydrodynamical models on stellar spectroscopy is evaluated by 
means of a differential 3D— ID abundance analysis based on a simple curve- of- growth 
method. The abundance of the trace element is varied independently in the 3D and ID 
calculations to match the equivalent width of a given spectral line; the difference between 
the two abundances represents the differential 3D— ID correction to the abundance of 
the trace element derived from that particular spectral feature. 

4. Results 

In this section, I will present some representative results of differential 3D— ID LTE 
abundance analyses at very low metallicity, where the impact of 3D models on the 
strengths of spectral lines is anticipated to be largest. Figure 3 shows the differential 
3D— ID LTE corrections to the Fe abundance for a very metal-poor red giant as derived 
from "fictitious" Fe I lines at A = 5000 A with varying equivalent widths (W\ between 
10 and 80 mA) and lower-level excitation potentials. § Such corrections are overall large 
and negative, particularly in the case of low-excitation lines, for which the differences 
between the predicted 3D and ID LTE Fe abundances are of the order of —1 dex. 
The behaviour of the differential 3D— ID LTE iron abundances from Fe I lines at low 
metallicity can be qualitatively explained by comparing the variations of the fraction of 
neutral-to-total Fe number densities with optical depth in the 3D hydrodynamical and in 
the corresponding ID marcs model atmospheres (Fig. 4). Under he assumption of LTE, 
in the ID MARCS model, iron is nearly completely ionized throughout the atmosphere; 
on the contrary, in the 3D model, the lower surface temperatures encountered in the 
upper photosphere allow a significant fraction of iron to recombine into neutral form. 
Hence, at a given Fe abundance, the density of absorbing neutral iron particles will be 
higher in the upper photosphere of the 3D model than in the ID model and Fe I lines 
will also appear stronger. Reversely, comparing with the ID case, a lower Fe abundance 
is required in 3D to match the same equivalent width of a given Fe I line. 

A similar behaviour holds for the differential 3D— ID LTE abundance corrections 

§ Fictitious lines (Stcffcn & Holweger 2002, Asplund 2005, Collet et al. 2006, Collet et al. 2007) provide 
a benchmark to analyse the behaviour of spectral lines solely as a function of lower-level excitation 
potential, wavelength, and strength, disregarding complications introduced by blends and wavelength 
dependency of continuous opacities in samples of real lines. 
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Figure 3. Differences between the Fe abundances derived in LTE from "fictitious" 
Fe I lines at A = 5000 A with a 3D hydrodynamical and a ID MARCS model atmosphere 
of very metal-poor red giant. The differential abundances are plotted as a function 
of equivalent width W\ and lower-level excitation potential of the lines, and for two 
different choices of the micro-turbulence parameter £ in the ID calculations. 



derived from lines of other neutral metals (Asplund 2005, Collet et al. 2007). Molecule 
formation as well shows an extreme sensitivity to temperature in the upper photospheric 
layers of late-type stars. This temperature sensitivity causes the 3D— ID LTE corrections 
to elemental abundances derived from the analysis of molecular features to also be 
large and negative in very metal-poor stars (Asplund & Garcia Perez 2001, Collet 
et al. 2006, Collet et al. 2007). 

5. 3D abundance analysis of extreme halo stars 

An interesting application of 3D model atmospheres of very metal-poor stars is 
the investigation of the impact of stellar granulation on the abundance analysis of 
the two extremely iron-poor halo stars HE0107-5240 and HE1327-2326 (Christlieb 
et al. 2002, Frebel et al. 2005). These two stars are remarkable in the sense that the 
abundance of the iron-peak elements in their atmospheres is the lowest ever observed 
in stellar environments ([Fe/H]< —5), while at the same time they are characterized by 
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Figure 4. Ratio of neutral to total iron number densities (riFei/^Fei) as a function 
of optical depth in the continuum at 5000 A in the atmosphere of a red giant star at 
[Fe/H]=— 3. Gray shaded area: distribution of riFei/^Fei values predicted with the use 
of a 3D model atmosphere: darker areas indicate values with higher occurrence. Over- 
plotted is the fraction of neutral iron as a function of optical depth in the corresponding 
ID MARCS model atmosphere (solid line). 



very large over-abundances of carbon, nitrogen, and oxygen with respect to iron. The 
interest aroused by these two objects comes from the consideration that HE0107— 5240 
and HE1327— 2326 might be direct descendants of a previous generation of metal-free 
stars. Various hypotheses have been proposed to explain the origin of the two stars but 
in order to identify the most plausible formation scenario, an accurate determination of 
the their chemical composition is necessary. Collet et al. (2006) and Frebel et al. (2008) 
have performed an abundance analysis of HE0 107— 5240 and HE1327— 2326 with the 
aid of 3D model atmospheres and following the basic procedure sketched in Sec. 3. The 
main results of the 3D — ID LTE abundance analysis of the two stars are summarized 
in Tab. 1. The 3D LTE Fe abundance is estimated to be about —0.2 dex lower than 
according to the ID analysis, while the abundances of carbon, nitrogen, and oxygen 
derived from CH, NH, and OH molecular lines are decreased by —0.7 dex or more. In 
the case of HE0107-5240, it is possible to derive the value of C and N abundances from 
more than one molecular indicator. Interestingly, while Christlieb et al. (2004) found 
a discrepancy of about 0.3 dex between the ID carbon abundance values derived from 



Three-dimensional models of metal-poor stars 



9 



Table 1. Average differential 3D— ID LTE corrections to the C, N, and O abundances 
in the extremely iron-poor halo stars HEOf 07— 5240 and HEf 327— 2326 as derived from 
low-excitation molecular lines; average corrections to the Fe abundance estimated from 
Fe I lines are also given. Note: abundances are expressed in the customary logarithmic 
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-2326 


logef D 
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lo g e 3D 
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CH 


6.81 


5.75 


6.90 


6.21 
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7.11 
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NH 


4.83 


3.8 


6.79 


6.10 
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CN 


5.22 (4.93) d 


3.2 









OH 


5.66 


5.0 


6.84 


6.12 


Fe 


Fe i 


2.06 


1.8 


1.79 


1.5 



a (Christlieb et al. 2004, Bessell et al. 2004, Bessell & Christlieb 2005) 

b (Collet et al. 2006) 

c (Frebel et al. 2008) 

d Assuming logei D (C) = 6.81 (7.11) 

CH and C2 lines, the 3D analysis brings the abundances from these two indicators down 
to the same value of loge(C)~5.7 dex. The 3D analysis of CN lines in HE0107-5240, 
however, returns an extremely large and negative 3D— ID correction to the nitrogen 
abundance (about —2 dex), which is not consistent with the 3D— ID correction evaluated 
from NH lines (see Tab. 1). This discrepancy could be, on the other hand, ascribable to 
the highly uncertain g /-values of NH lines (Collet et al. 2006). This problem certainly 
deserves further investigation. 

6. Discussion 

The prediction of photospheric temperatures significantly below the radiative 
equilibrium value is a central result of surface convection simulations of metal-poor late- 
type stars. One implication of this result is that 3D— ID abundance corrections in very 
metal-poor stars are in general significantly larger than other systematic errors usually 
quoted in classical abundance analyses (Gustafsson 2004). It is therefore necessary to 
investigate to what extent the cooler temperature stratification of the 3D models is 
controlled by the assumptions made during the construction of the simulations. 

A relevant question in this respect is whether the opacity-binning scheme used 
for the calculation of the radiative heating rates represents a reasonably accurate 
approximation to the solution of the radiative transfer equation for the full set of 
wavelengths of the original opacity distribution functions (ODFs). Figure 5 shows the 
results of test calculations where the radiative heating rates are computed with the two 
different approaches in "1.5D" approximation along individual columns in a vertical slice 
of a metal-poor red giant simulation snapshot. The correlation between the opacity- 
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Figure 5. Comparison between the radiative heating rates computed in "1.5D" 
approximation for all grid-points on a vertical slice of a simulation snapshot of very 
metal-poor red giant using the opacity-binning scheme ( "BINS" ) and monochromatic 
radiative transfer with opacity distribution functions ( "ODFs" ) . The radiative heating 
rates per unit mass (Q ra d) are normalized with respect to the specie heat (per unit 
mass) at constant pressure (Cp ). A magnied view of the plot in the region of low 
radiative heating rates is shown in the small panel. 



binning and ODF-based schemes is very good, suggesting that the former is probably 
accurate enough in reproducing radiative heating and temperatures at the surface of 
metal-poor 3D simulations. 

Another important issue is whether the treatment of scattering as true absorption 
in the solution of the radiative transfer equation introduces systematic errors in the 
temperature at the surface of 3D convection simulations. According to preliminary 
tests with ID OS-marcs (Gustafsson et al. 2008) model atmospheres of metal-poor 
giants, the ID temperature stratifications predicted when scattering is included as true 
absorption is overall hotter than in calculations where scattering is treated correctly 
as such. In the example considered in Fig. 6, temperature differences between the two 
cases reach 300 K in the upper photosphere. This result would therefore suggest that 
the radiative heating rates and temperature might actually be overestimated at the 
surface of 3D models of metal-poor stars. It is not obvious, however, that the effect 
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Figure 6. Comparison between the temperature stratifications of two OS-MARCS 
model metal-poor giant atmospheres computed for the same stellar parameters but 
with two different implementations of continuous scattering. Solid line: scattering 
correctly treated as such in the source function (standard version of OS-MARCS) ; Dashed 
line: scattering included as true absorption, as in the 3D simulations presented here. 



of implementing scattering as true absorption would proceed in the same direction in 
3D. A full 3D treatment of scattering is paramount in this respect to establish whether 
the surface temperatures possibly underestimated in present metal-poor 3D models. I 
refer to the contribution of W. Hayek in these proceedings for more details about 3D 
radiative transfer with scattering in convection simulations. 

It is also crucial to assess how spectral line formation is affected by the 
approximations involved in the solution of the radiative transfer equation. The 
treatment of scattering as true absorption is a matter of concern also for 3D line 
formation calculations. Rayleigh scattering off H I is, in fact, an important source 
of extinction in the UV and blue part of the spectrum. At those wavelength, 
implementing scattering as true absorption causes the outgoing flux in the continuum 
to be underestimated and hence leads to the prediction of weaker spectral lines. The 
effect is expected to be more important at very low metallicities, due to the weak line- 
blocking, and in in metal-poor 3D model atmospheres where the density of scatterers 
(H i) is high because of the low surface temperatures. The use of a differential 3D— ID 
abundance analysis ensures at least that the uncertainties in the treatment of scattering 
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Table 2. Non-LTE Fe abundances derived for the extreme halo stars HE0107-5240 
and HE1327— 2326 using ID MARCS models and the mean temperature stratification 
versus optical depth from 3D models. The non-LTE calculations assume the model 
Fe atom by Collet et al. (2005), efficient Drawin-like H+Fe collisions (Drawin 1968, 
Drawin 1969), and thermalization of the uppermost Fe I levels. 





HE0107-5240 


HE1327-2326 


Model 


[Fe/H] LTE [Fe/H] non _LTE 


[Fe/H] L TE [Fe/H] non _LTE 


ID MARCS 


-5.40 -4.65 


-5.70 -5.15 


mean 3D 


-5.60 -4.75 


-5.95 -5.05 



are minimized. 

Finally, I would like to caution that many of the abundance indicators considered 
in the present contribution (e.g. neutral metals and possibly molecules as well) are 
most likely affected by departures from LTE or from chemical equilibrium at the local 
temperature. Neutral iron, for instance, is expected to be prone to such departures. 
The main non-LTE mechanism for Fe I in the photospheres of late-type stars is namely 
efficient over-ionization driven by the UV radiation field coming from layers deeper in. 
This causes Fe I levels to be underpopulated with respect to LTE, makes Fe I lines weaker 
than the LTE case at a given Fe abundance, and, consequently, leads to the prediction 
of higher Fe abundances. Fe I departures from LTE have been estimated by (Collet 
et al. 2006) for the two extreme halo stars HE0107-5240 and HE1327-2326 by means 
of a ID analysis (Tab. 2). The results of this preliminary calculations clearly suggest 
that the non-LTE effects are considerable and opposite to the corrections for stellar 
granulation. This implies that a combined 3D non-LTE analysis of Fe I line formation 
in a similar way as done by Asplund et al. (2003) and Barklem et al. (2003) for Li I is 
essential to produce realistic Fe abundance determinations. For completeness, I should 
also add that the departures of Fe I from LTE depend crucially on the magnitude of 
H+Fe collisional excitation and ionization processes, whose cross-sections are still poorly 
known. This explains, at least in part, why current estimates of non-LTE effects on Fe I 
line formation are still hampered by large uncertainties (see also review by A. Korn in 
these proceedings). 

7. Conclusions 

I have presented here some illustrative results of the application of 3D surface 
convection simulations of very metal-poor stars to spectral line formation in LTE. The 
differences between the predicted temperature stratifications of 3D hydrodynamical 
simulations and ID hydrostatic model atmospheres, and the 3D temperature and 
density inhomogeneities and correlated velocity gradients can significantly affect line 
strengths and, in turn, the elemental abundances inferred from spectral lines. At very 
low metallicities, the deviations of the mean 3D thermal structure from the classical 
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ID stratification are largest, and cause the 3D— ID LTE differential abundances to 
be negative and considerable for lines of neutral species (down to about —1 dex for 
weak low-excitation Fe I line). Corrections to CNO abundances derived from weak 
low-excitation CH, NH, and OH features are also found to be typically in the range 
—0.5 dex to —1.0 dex in very metal-poor giants (Collet et al. 2007). Finally I have 
discussed possible systematic errors affecting the present 3D abundance analyses, such 
as departures of Fe I line formation from LTE; non-LTE corrections to the Fe abundance 
are opposite to and, according to preliminary ID test calculations, of the same order of 
magnitude as the ones due to stellar granulation. 
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